A Bayesian framework to study tumor subclone-specific expression by combining bulk DNA and single-cell RNA sequencing data.

Genetic and gene expression heterogeneity is an essential hallmark of many tumors, allowing the cancer to evolve and to develop resistance to treatment. Currently, the most commonly used data types for studying such heterogeneity are bulk tumor/normal whole-genome or whole-exome sequencing (WGS, WES); and single-cell RNA sequencing (scRNA-seq), respectively. However, tools are currently lacking to link genomic tumor subclonality with transcriptomic heterogeneity by integrating genomic and single-cell transcriptomic data collected from the same tumor. To address this gap, we developed scBayes, a Bayesian probabilistic framework that uses tumor subclonal structure inferred from bulk DNA sequencing data to determine the subclonal identity of cells from single-cell gene expression (scRNA-seq) measurements. Grouping together cells representing the same genetically defined tumor subclones allows comparison of gene expression across different subclones, or investigation of gene expression changes within the same subclone across time (i.e., progression, treatment response, or relapse) or space (i.e., at multiple metastatic sites and organs). We used simulated data sets, in silico synthetic data sets, as well as biological data sets generated from cancer samples to extensively characterize and validate the performance of our method, as well as to show improvements over existing methods. We show the validity and utility of our approach by applying it to published data sets and recapitulating the findings, as well as arriving at novel insights into cancer subclonal expression behavior in our own data sets. We further show that our method is applicable to a wide range of single-cell sequencing technologies including single-cell DNA sequencing as well as Smart-seq and 10x Genomics scRNA-seq protocols.

Inhibition of TRPA1, Endoplasmic Reticulum Stress, Human Airway Epithelial Cell Damage, and Ectopic MUC5AC Expression by Vasaka (Adhatoda vasica; Malabar Nut) Tea.

This study tested whether a medicinal plant, Vasaka, typically consumed as a tea to treat respiratory malaise, could protect airway epithelial cells (AECs) from wood smoke particle-induced damage and prevent pathological mucus production. Wood/biomass smoke is a pneumotoxic air pollutant. Mucus normally protects the airways, but excessive production can obstruct airflow and cause respiratory distress. Vasaka tea pre- and co-treatment dose-dependently inhibited mucin 5AC (MUC5AC) mRNA induction by AECs treated with wood smoke particles. This correlated with transient receptor potential ankyrin-1 (TRPA1) inhibition, an attenuation of endoplasmic reticulum (ER) stress, and AEC damage/death. Induction of mRNA for anterior gradient 2, an ER chaperone/disulfide isomerase required for MUC5AC production, and TRP vanilloid-3, a gene that suppresses ER stress and wood smoke particle-induced cell death, was also attenuated. Variable inhibition of TRPA1, ER stress, and MUC5AC mRNA induction was observed using selected chemicals identified in Vasaka tea including vasicine, vasicinone, apigenin, vitexin, isovitexin, isoorientin, 9-oxoODE, and 9,10-EpOME. Apigenin and 9,10-EpOME were the most cytoprotective and mucosuppressive. Cytochrome P450 1A1 (CYP1A1) mRNA was also induced by Vasaka tea and wood smoke particles. Inhibition of CYP1A1 enhanced ER stress and MUC5AC mRNA expression, suggesting a possible role in producing protective oxylipins in stressed cells. The results provide mechanistic insights and support for the purported benefits of Vasaka tea in treating lung inflammatory conditions, raising the possibility of further development as a preventative and/or restorative therapy.

MITF Is Regulated by Redox Signals Controlled by the Selenoprotein Thioredoxin Reductase 1.

TR1 and other selenoproteins have paradoxical effects in melanocytes and melanomas. Increasing selenoprotein activity with supplemental selenium in a mouse model of UV-induced melanoma prevents oxidative damage to melanocytes and delays melanoma tumor formation. However, TR1 itself is positively associated with progression in human melanomas and facilitates metastasis in melanoma xenografts. Here, we report that melanocytes expressing a microRNA directed against TR1 (TR1low) grow more slowly than control cell lines and contain significantly less melanin. This phenotype is associated with lower tyrosinase (TYR) activity and reduced transcription of tyrosinase-like protein-1 (TYRP1). Melanoma cells in which the TR1 gene (TXNRD1) was disrupted using Crispr/Cas9 showed more dramatic effects including the complete loss of the melanocyte-specific isoform of MITF; other MITF isoforms were unaffected. We provide evidence that TR1 depletion results in oxidation of MITF itself. This newly discovered mechanism for redox modification of MITF has profound implications for controlling both pigmentation and tumorigenesis in cells of the melanocyte lineage.

Thioredoxin Reductase 1 Modulates Pigmentation and Photobiology of Murine Melanocytes in vivo.

Pigment-producing melanocytes overcome frequent oxidative stress in their physiological role of protecting the skin against the deleterious effects of solar UV irradiation. This is accomplished by the activity of several endogenous antioxidant systems, including the thioredoxin antioxidant system, in which thioredoxin reductase 1 (TR1) plays an important part. To determine whether TR1 contributes to the redox regulation of melanocyte homeostasis, we have generated a selective melanocytic Txnrd1-knockout mouse model (Txnrd1mel‒/‒), which exhibits a depigmentation phenotype consisting of variable amelanotic ventral spotting and reduced pigmentation on the extremities (tail tip, ears, and paws). The antioxidant role of TR1 was further probed in the presence of acute neonatal UVB irradiation, which stimulates melanocyte activation and introduces a spike in oxidative stress in the skin microenvironment. Interestingly, we observed a significant reduction in overall melanocyte count and proliferation in the absence of TR1. Furthermore, melanocytes exhibited an elevated level of UV-induced DNA damage in the form of 8-oxo-2'-deoxyguanosine after acute UVB treatment. We also saw an engagement of compensatory antioxidant mechanisms through increased nuclear localization of transcription factor NRF2. Altogether, these data indicate that melanocytic TR1 positively regulates melanocyte homeostasis and pigmentation during development and protects against UVB-induced DNA damage and oxidative stress.

Topsentinol L Trisulfate, a Marine Natural Product That Targets Basal-like and Claudin-Low Breast Cancers.

Patients diagnosed with basal-like breast cancer suffer from poor prognosis and limited treatment options. There is an urgent need to identify new targets that can benefit patients with basal-like and claudin-low (BL-CL) breast cancers. We screened fractions from our Marine Invertebrate Compound Library (MICL) to identify compounds that specifically target BL-CL breast cancers. We identified a previously unreported trisulfated sterol, i.e., topsentinol L trisulfate (TLT), which exhibited increased efficacy against BL-CL breast cancers relative to luminal/HER2+ breast cancer. Biochemical investigation of the effects of TLT on BL-CL cell lines revealed its ability to inhibit activation of AMP-activated protein kinase (AMPK) and checkpoint kinase 1 (CHK1) and to promote activation of p38. The importance of targeting AMPK and CHK1 in BL-CL cell lines was validated by treating a panel of breast cancer cell lines with known small molecule inhibitors of AMPK (dorsomorphin) and CHK1 (Ly2603618) and recording the increased effectiveness against BL-CL breast cancers as compared with luminal/HER2+ breast cancer. Finally, we generated a drug response gene-expression signature and projected it against a human tumor panel of 12 different cancer types to identify other cancer types sensitive to the compound. The TLT sensitivity gene-expression signature identified breast and bladder cancer as the most sensitive to TLT, while glioblastoma multiforme was the least sensitive.

A `one-two punch' therapy strategy to target chemoresistance in estrogen receptor positive breast cancer.

Cancer cell phenotypes evolve during a tumor's treatment. In some cases, tumor cells acquire cancer stem cell-like (CSL) traits such as resistance to chemotherapy and diminished differentiation; therefore, targeting these cells may be therapeutically beneficial. In this study we show that in progressive estrogen receptor positive (ER+) metastatic breast cancer tumors, resistant subclones that emerge following chemotherapy have increased CSL abundance. Further, in vitro organoid growth of ER+ patient cancer cells also shows that chemotherapy treatment leads to increased abundance of ALDH+/CD44+ CSL cells. Chemotherapy induced CSL abundance is blocked by treatment with a pan-HDAC inhibitor, belinostat. Belinostat treatment diminished both mammosphere formation and size following chemotherapy, indicating a decrease in progenitor CSL traits. HDAC inhibitors specific to class IIa (HDAC4, HDAC5) and IIb (HDAC6) were shown to primarily reverse the chemo-resistant CSL state. Single-cell RNA sequencing analysis with patient samples showed that HDAC targets and MYC signaling were promoted by chemotherapy and inhibited upon HDAC inhibitor treatment. In summary, HDAC inhibition can block chemotherapy-induced drug resistant phenotypes with 'one-two punch' strategy in refractory breast cancer cells.

Exploiting collateral sensitivity controls growth of mixed culture of sensitive and resistant cells and decreases selection for resistant cells in a cell line model.

BACKGROUND: CDK4/6 inhibitors such as ribociclib are becoming widely used targeted therapies in hormone-receptor-positive (HR+) human epidermal growth factor receptor 2-negative (HER2-) breast cancer. However, cancers can advance due to drug resistance, a problem in which tumor heterogeneity and evolution are key features. METHODS: Ribociclib-resistant HR+/HER2- CAMA-1 breast cancer cells were generated through long-term ribociclib treatment. Characterization of sensitive and resistant cells were performed using RNA sequencing and whole exome sequencing. Lentiviral labeling with different fluorescent proteins enabled us to track the proliferation of sensitive and resistant cells under different treatments in a heterogeneous, 3D spheroid coculture system using imaging microscopy and flow cytometry. RESULTS: Transcriptional profiling of sensitive and resistant cells revealed the downregulation of the G2/M checkpoint in the resistant cells. Exploiting this acquired vulnerability; resistant cells exhibited collateral sensitivity for the Wee-1 inhibitor, adavosertib (AZD1775). The combination of ribociclib and adavosertib achieved additional antiproliferative effect exclusively in the cocultures compared to monocultures, while decreasing the selection for resistant cells. CONCLUSIONS: Our results suggest that optimal antiproliferative effects in heterogeneous cancers can be achieved via an integrative therapeutic approach targeting sensitive and resistant cancer cell populations within a tumor, respectively.

Circulating immune cell phenotype dynamics reflect the strength of tumor-immune cell interactions in patients during immunotherapy.

The extent to which immune cell phenotypes in the peripheral blood reflect within-tumor immune activity prior to and early in cancer therapy is unclear. To address this question, we studied the population dynamics of tumor and immune cells, and immune phenotypic changes, using clinical tumor and immune cell measurements and single-cell genomic analyses. These samples were serially obtained from a cohort of advanced gastrointestinal cancer patients enrolled in a trial with chemotherapy and immunotherapy. Using an ecological population model, fitted to clinical tumor burden and immune cell abundance data from each patient, we find evidence of a strong tumor-circulating immune cell interaction in responder patients but not in those patients that progress on treatment. Upon initiation of therapy, immune cell abundance increased rapidly in responsive patients, and once the peak level is reached tumor burden decreases, similar to models of predator-prey interactions; these dynamic patterns were absent in nonresponder patients. To interrogate phenotype dynamics of circulating immune cells, we performed single-cell RNA sequencing at serial time points during treatment. These data show that peripheral immune cell phenotypes were linked to the increased strength of patients' tumor-immune cell interaction, including increased cytotoxic differentiation and strong activation of interferon signaling in peripheral T cells in responder patients. Joint modeling of clinical and genomic data highlights the interactions between tumor and immune cell populations and reveals how variation in patient responsiveness can be explained by differences in peripheral immune cell signaling and differentiation soon after the initiation of immunotherapy.

Time- and dose-dependent gene expression analysis of macrophage response as a function of porosity of silica nanoparticles.

There is a limited amount of information available on gene expression regulation of macrophages in response to changing the time of exposure, concentration, and physicochemical properties of nanomaterials. In this study, RAW264.7 macrophages were treated with spherical nonporous and mesoporous silica nanoparticles of similar size at different incubation times and concentrations. RNA-sequencing was used to study transcriptional profiles. Bioinformatics analyses, functional annotation clustering, and network analyses were employed to understand signaling pathways of cellular response as a function of porosity, incubation time, and concentration. Porosity introduced drastic changes to the genomic response of macrophages at equitoxic concentrations and incubation times. Direct relations between increases in time and concentration with an increased number of differentially expressed genes were observed.

Genotoxicity of amorphous silica nanoparticles: Status and prospects.

Amorphous silica nanoparticles (SNPs) are widely used in biomedical applications and consumer products. Little is known, however, about their genotoxicity and potential to induce gene expression regulation. Despite recent efforts to study the underlying mechanisms of genotoxicity of SNPs, inconsistent results create a challenge. A variety of factors determine particle-cell interactions and underlying mechanisms. Further, high-throughput studies are required to carefully assess the impact of silica nanoparticle physicochemical properties on induction of genotoxic response in different cell lines and animal models. In this article, we review the strategies available for evaluation of genotoxicity of nanoparticles (NPs), survey current status of silica nanoparticle gene alteration and genotoxicity, discuss particle-mediated inflammation as a contributing factor to genotoxicity, identify existing gaps and suggest future directions for this research.

Differential gene expression patterns in vein regions susceptible versus resistant to neointimal hyperplasia.

Arteriovenous hemodialysis graft (AVG) stenosis results in thrombosis and AVG failure, but prevention of stenosis has been unsuccessful due in large part to our limited understanding of the molecular processes involved in neointimal hyperplasia (NH) formation. AVG stenosis develops chiefly as a consequence of highly localized NH formation in the vein-graft anastomosis region. Surprisingly, the vein region just downstream of the vein-graft anastomosis (herein termed proximal vein region) is relatively resistant to NH. We hypothesized that the gene expression profiles of the NH-prone and NH-resistant regions will be different from each other after graft placement, and analysis of their genomic profiles may yield potential therapeutic targets to prevent AVG stenosis. To test this, we evaluated the vein-graft anastomosis (NH-prone) and proximal vein (NH-resistant) regions in a porcine model of AVG stenosis with a porcine microarray. Gene expression changes in these two distinct vein regions, relative to the gene expression in unoperated control veins, were examined at early (5 days) and later (14 days) time points following graft placement. Global genomic changes were much greater in the NH-prone region than in the NH-resistant region at both time points. In the NH-prone region, genes related to regulation of cell proliferation and osteo-/chondrogenic vascular remodeling were most enriched among the significantly upregulated genes, and genes related to smooth muscle phenotype were significantly downregulated. These results provide insights into the spatial and temporal genomic modulation underlying NH formation in AVG and suggest potential therapeutic strategies to prevent and/or limit AVG stenosis.

Publisher Correction: Combating subclonal evolution of resistant cancer phenotypes.

The originally published version of this Article contained an error in Figure 4. In panel a, grey boxes surrounding the subclones associated with patients #2 and #4 obscured adjacent portions of the heatmap. This error has now been corrected in both the PDF and HTML versions of the Article.

Global gene expression analysis of macrophage response induced by nonporous and porous silica nanoparticles.

Little is known about the global gene expression profile of macrophages in response to changes in size and porosity of silica nanoparticles (SNPs). Spherical nonporous SNPs of two different diameters, and mesoporous spherical SNPs with comparable size were characterized. Reactive oxygen species, mitochondrial membrane potential, lysosome degradation capacity, and lysosome pH were measured to evaluate the influence of nonporous and mesoporous SNPs on mitochondrial and lysosomal function. RNA-sequencing was utilized to generate transcriptional profiles of RAW264.7 macrophages exposed to non-toxic SNP doses. DESeq2, limma, and BinReg2 software were used to analyze the data based on both unsupervised and supervised strategies to identify genes with greatest differences among NP treatments. Utilizing GATHER and DAVID software, possible induced pathways were studied. We found that mesoporous silica nanoparticles are capable of altering gene expression in macrophages at doses that do not elicit acute cytotoxicity, while gene transcription was minimally affected by nonporous SNPs.

Combating subclonal evolution of resistant cancer phenotypes.

Metastatic breast cancer remains challenging to treat, and most patients ultimately progress on therapy. This acquired drug resistance is largely due to drug-refractory sub-populations (subclones) within heterogeneous tumors. Here, we track the genetic and phenotypic subclonal evolution of four breast cancers through years of treatment to better understand how breast cancers become drug-resistant. Recurrently appearing post-chemotherapy mutations are rare. However, bulk and single-cell RNA sequencing reveal acquisition of malignant phenotypes after treatment, including enhanced mesenchymal and growth factor signaling, which may promote drug resistance, and decreased antigen presentation and TNF-α signaling, which may enable immune system avoidance. Some of these phenotypes pre-exist in pre-treatment subclones that become dominant after chemotherapy, indicating selection for resistance phenotypes. Post-chemotherapy cancer cells are effectively treated with drugs targeting acquired phenotypes. These findings highlight cancer's ability to evolve phenotypically and suggest a phenotype-targeted treatment strategy that adapts to cancer as it evolves.

Integrative analyses reveal signaling pathways underlying familial breast cancer susceptibility.

The signaling events that drive familial breast cancer (FBC) risk remain poorly understood. While the majority of genomic studies have focused on genetic risk variants, known risk variants account for at most 30% of FBC cases. Considering that multiple genes may influence FBC risk, we hypothesized that a pathway-based strategy examining different data types from multiple tissues could elucidate the biological basis for FBC. In this study, we performed integrated analyses of gene expression and exome-sequencing data from peripheral blood mononuclear cells and showed that cell adhesion pathways are significantly and consistently dysregulated in women who develop FBC. The dysregulation of cell adhesion pathways in high-risk women was also identified by pathway-based profiling applied to normal breast tissue data from two independent cohorts. The results of our genomic analyses were validated in normal primary mammary epithelial cells from high-risk and control women, using cell-based functional assays, drug-response assays, fluorescence microscopy, and Western blotting assays. Both genomic and cell-based experiments indicate that cell-cell and cell-extracellular matrix adhesion processes seem to be disrupted in non-malignant cells of women at high risk for FBC and suggest a potential role for these processes in FBC development.

Gene-expression patterns in peripheral blood classify familial breast cancer susceptibility.

BACKGROUND: Women with a family history of breast cancer face considerable uncertainty about whether to pursue standard screening, intensive screening, or prophylactic surgery. Accurate and individualized risk-estimation approaches may help these women make more informed decisions. Although highly penetrant genetic variants have been associated with familial breast cancer (FBC) risk, many individuals do not carry these variants, and many carriers never develop breast cancer. Common risk variants have a relatively modest effect on risk and show limited potential for predicting FBC development. As an alternative, we hypothesized that additional genomic data types, such as gene-expression levels, which can reflect genetic and epigenetic variation, could contribute to classifying a person's risk status. Specifically, we aimed to identify common patterns in gene-expression levels across individuals who develop FBC. METHODS: We profiled peripheral blood mononuclear cells from women with a family history of breast cancer (with or without a germline BRCA1/2 variant) and from controls. We used the support vector machines algorithm to differentiate between patients who developed FBC and those who did not. Our study used two independent datasets, a training set of 124 women from Utah (USA) and an external validation (test) set from Ontario (Canada) of 73 women (197 total). We controlled for expression variation associated with clinical, demographic, and treatment variables as well as lymphocyte markers. RESULTS: Our multigene biomarker provided accurate, individual-level estimates of FBC occurrence for the Utah cohort (AUC = 0.76 [0.67-84]) . Even at their lower confidence bounds, these accuracy estimates meet or exceed estimates from alternative approaches. Our Ontario cohort resulted in similarly high levels of accuracy (AUC = 0.73 [0.59-0.86]), thus providing external validation of our findings. Individuals deemed to have "high" risk by our model would have an estimated 2.4 times greater odds of developing familial breast cancer than individuals deemed to have "low" risk. CONCLUSIONS: Together, these findings suggest that gene-expression levels in peripheral blood cells reflect genomic variation associated with breast cancer risk and that such data have potential to be used as a non-invasive biomarker for familial breast cancer risk.

The role of thioredoxin reductase 1 in melanoma metabolism and metastasis.

Although significant progress has been made in targeted and immunologic therapeutics for melanoma, many tumors fail to respond, and most eventually progress when treated with the most efficacious targeted combination therapies thus far identified. Therefore, alternative approaches that exploit distinct melanoma phenotypes are necessary to develop new approaches for therapeutic intervention. Tissue microarrays containing human nevi and melanomas were used to evaluate levels of the antioxidant protein thioredoxin reductase 1 (TR1), which was found to increase as a function of disease progression. Melanoma cell lines revealed metabolic differences that correlated with TR1 levels. We used this new insight to design a model treatment strategy that creates a synthetic lethal interaction wherein targeting TR1 sensitizes melanoma to inhibition of glycolytic metabolism, resulting in a decrease in metastases in vivo. This approach holds the promise of a new clinical therapeutic strategy, distinct from oncoprotein inhibition.

Genomic classification of the RAS network identifies a personalized treatment strategy for lung cancer.

Better approaches are needed to evaluate a single patient's drug response at the genomic level. Targeted therapy for signaling pathways in cancer has met limited success in part due to the exceedingly interwoven nature of the pathways. In particular, the highly complex RAS network has been challenging to target. Effectively targeting the pathway requires development of techniques that measure global network activity to account for pathway complexity. For this purpose, we used a gene-expression-based biomarker for RAS network activity in non-small cell lung cancer (NSCLC) cells, and screened for drugs whose efficacy was significantly highly correlated to RAS network activity. Results identified EGFR and MEK co-inhibition as the most effective treatment for RAS-active NSCLC amongst a panel of over 360 compounds and fractions. RAS activity was identified in both RAS-mutant and wild-type lines, indicating broad characterization of RAS signaling inclusive of multiple mechanisms of RAS activity, and not solely based on mutation status. Mechanistic studies demonstrated that co-inhibition of EGFR and MEK induced apoptosis and blocked both EGFR-RAS-RAF-MEK-ERK and EGFR-PI3K-AKT-RPS6 nodes simultaneously in RAS-active, but not RAS-inactive NSCLC. These results provide a comprehensive strategy to personalize treatment of NSCLC based on RAS network dysregulation and provide proof-of-concept of a genomic approach to classify and target complex signaling networks.

A chimeric p53 evades mutant p53 transdominant inhibition in cancer cells.

Because of the dominant negative effect of mutant p53, there has been limited success with wild-type (wt) p53 cancer gene therapy. Therefore, an alternative oligomerization domain for p53 was investigated to enhance the utility of p53 for gene therapy. The tetramerization domain of p53 was substituted with the coiled-coil (CC) domain from Bcr (breakpoint cluster region). Our p53 variant (p53-CC) maintains proper nuclear localization in breast cancer cells detected via fluorescence microscopy and shows a similar expression profile of p53 target genes as wt-p53. Additionally, similar tumor suppressor activities of p53-CC and wt-p53 were detected by terminal deoxynucleotidyl transferase dUTP nick end labeling (TUNEL), annexin-V, 7-aminoactinomycin D (7-AAD), and colony-forming assays. Furthermore, p53-CC was found to cause apoptosis in four different cancer cell lines, regardless of endogenous p53 status. Interestingly, the transcriptional activity of p53-CC was higher than wt-p53 in 3 different reporter gene assays. We hypothesized that the higher transcriptional activity of p53-CC over wt-p53 was due to the sequestration of wt-p53 by endogenous mutant p53 found in cancer cells. Co-immunoprecipitation revealed that wt-p53 does indeed interact with endogenous mutant p53 via its tetramerization domain, while p53-CC escapes this interaction. Therefore, we investigated the impact of the presence of a transdominant mutant p53 on tumor suppressor activities of wt-p53 and p53-CC. Overexpression of a potent mutant p53 along with wt-p53 or p53-CC revealed that, unlike wt-p53, p53-CC retains the same level of tumor suppressor activity. Finally, viral transduction of wt-p53 and p53-CC into a breast cancer cell line that harbors a tumor derived transdominant mutant p53 validated that p53-CC indeed evades sequestration and consequent transdominant inhibition by endogenous mutant p53.

Transcriptional responses of human aortic endothelial cells to nanoconstructs used in biomedical applications.

Understanding the potential toxicities of manufactured nanoconstructs used for drug delivery and biomedical applications may help improve their safety. We sought to determine if surface-modified silica nanoparticles and poly(amido amine) dendrimers elicit genotoxic responses on vascular endothelial cells. The nanoconstructs utilized in this study had a distinct geometry (spheres vs worms) and surface charge, which were used to evaluate the contributions of these parameters to any potential adverse effects of these materials. Time-dependent cytotoxicity was found for surfaced-functionalized but geometrically distinct silica materials, while amine-terminated dendrimers displayed time-independent cytotoxicity and carboxylated dendrimers were nontoxic in our assays. Transcriptomic evaluation of human aortic endothelial cell (HAEC) responses indicated time-dependent gene induction following silica exposure, consisting of cell cycle gene repression and pro-inflammatory gene induction. However, the dendrimers did not induce genomic toxicity, despite displaying general cytotoxicity.